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Abstract 

A systematic analysis of a continuous version of a binomial lattice, 
containing a real parameter 7 and covering the Toda field equation as 
7 00, is carried out in the framework of group theory. The symme- 
try algebra of the equation is derived. Reductions by one-dimensional 
and two-dimensional subalgebras of the symmetry algebra and their cor- 
responding subgroups, yield notable field equations in lower dimensions 
whose solutions allow to find exact solutions to the original equation. 
Some reduced equations turn out to be related to potentials of physical in- 
terest, such as the Fermi-Pasta-Ulam and the Killingbeck potentials, and 
others. An instanton-like approximate solution is also obtained which re- 
produces the Eguchi-Hanson instanton configuration for 7 — > 00. Further- 
more, the equation under consideration is extended to (n-|- l)-dimensions. 
A spherically symmetric form of this equation, studied by means of the 
symmetry approach, provides conformally invariant classes of field equa- 
tions comprising remarkable special cases. One of these (n = 4) enables 
us to establish a connection with the Euclidean Yang-Mills equations, an- 
other appears in the context of Differential Geometry in relation to the 
socalled Yamabe problem. All the properties of the reduced equations are 
shared by the spherically symmetric generalized field equation. 



1 



1 Introduction 



The continuous (or long-wave) approximations of chains of particles can provide 
a fruitful theoretical framework to be used as a guide in the study of the cor- 
responding original discrete systems. These models are generally formulated as 
nonlinear field equations. When such equations are integrable, exact solutions 
can be found as, for example, solitons, vortices and shock-wave. 

Continuous limits of lattice systems may be considered as phenomenological 
models with a proper identity, and their study could be interesting apart from 
the reduction procedure applied in passing to the continuous representation 
[0. This feature is shared by different models, among which the well-known 
continuous form 

Uxx + Uvv-k{e")zz^O (1.1) 

of the two-dimensional discrete Toda equation, where u ~ u{x, y, z), plays a 
basic role (Subscripts denote partial derivatives and k — ±1). 

Indeed, Eq. (1.1) appears in a variety of physical areas, running from the 
theory of Hamiltonian systems to general relativity |^ , . In the latter context, 
Eq. (1.1) occurs in the theory of self-dual Einstein spaces, where is rechristened 
heavenly equation Q , (|] . This emerges as a limit case from the Toda molecule 
equation which is exactly integrable 

Another nonlinear field equation associated with a nonlinear lattice which 
deserves to be investigated is 



7-1 



0, 



(1.2) 



where 7 is a (real) parameter such that 7^ 0, 1. 

Really, after suitable rescalings, Eq. (O) can be interpreted as the contin- 
uous limit of a uniform two-dimensional nonlinear lattice of N particles inter- 
acting through the nearest-neighbor potential 
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7 



(1 + 6„r„) 



(1.3) 



where a„ and 6„ are constants of the n-th nonlinear spring, r„ — Un — Vn 
and Tjn is the displacement of the n-particle from its equilibrium position. 



We shall call "binomial latti ce" the chain described by the potential (1.3). 

We notice that the function ( |l.3D covers the Toda potential (7 ^ cx)) ||^, the 
harmonic potential (7 = 2), and the potential used by Fermi, Pasta and Ulam 
in their computer experiment in the early's 1950 (7 = 3) |9| . 

As one expects, for 7 — > cx3 Eq. (1.2) becomes Eq. ( |l . 1| ) . 

The purpose of this paper is both to find the symmetry structure and exact 
solutions of Eq. (1.2). An effective method to get insight into the symmetry 
properties of Eq. (1.2) and, consequently, to construct explicit configurations, 
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is based on the reduction approach [Q, which exploits group-theoretical tech- 
niques. Applying this procedure, we obtain the symmetry algebra of Eq. (1.2). 
The generator of the corresponding symmetry group has been determined in 
part with the help of a computer, by means of the symbolic language MAPLE 
V Release 4 |11 . Our result is that Eq. ( |1.2| ) admits a /inzie-dimensional Lie 
group of symmetries, i.e. a local group G of transformations acting on the in- 
dependent variables {x, y, z) and the dependent variable u with the property 
that whenever u{x, y, z) is a solution of Eq. (L2), then u' ^ {g o u){x' , y', z') 
is also a solution for any g & G. 

We remind the reader that in contrast with what happens for Eq. (L2), the 
Toda equation handled within the group theory allows an infinite- 

dimcnsional symmetry algebra, a realization of which is given by generators 
obeying a Virasoro algebra without a central charge (Witt algebra). Further- 
more, certain reduced equations give rise to instanton and meron-like solutions 
endowed with integer and fractional topological numbers, respectively. 

Thus, the comparative analysis of the properties of Eqs. (LI) and ( |l.2| ) is of 
particular relevance, keeping in mind also the fact that both the equations come 
from the continuous approximation of physically significant lattice systems. 

Another important characteristic common to Eqs. (|l|) and (O ) is that 
they are related to nonlinear systems of the hydrodynamical type ^|. Although 
interesting, here this aspect will not be treated. 

The main results achieved in this article are presented in the Sections in 
which the paper is organized. 

In Sec. 2 we find the symmetry algebra of Eq. (1.2). It turns out that the 
related independent infinitesimal operators obey a closed set of commutation 
rules. Starting from special linear combinations of these operators, in Sec. 3 
the corresponding symmetry group transformations are derived. These provide 
reduced differential equations leading to exact solutions to Eq. ( LS ) . Notable 
cases arise for 7 — 3, —2, |, |. This choice of values is motivated by the fact 
that, in correspondence, the function (|l.3| ) coincides with potentials of physical 
interest, such as the Fermi-Pasta-Ulam (7 = 3) |^ and a potential involved in 
the Thomas-Fermi model of an atom (7 = 5/2) p. 116). We point out 

that for 7 = 3, a remarkable representation of the inverse Weierstrass function 
is given in terms of the Gauss hypergeometric function. In Sec. 4, we obtain 
an approximate instanton-like configuration, holding for large values of the 
parameter 7, which reproduces the Eguchi-Hanson instanton solution of Eq. 
( PI ) in the limit 7 ^ 00 |T1 . 



In Sec. 5 we extend Eq. (1.2) to the (n + l)-dimensional case. The explo- 
ration of the symmetry properties of this more general nonlinear field equation, 
that is Eq. (5.1), has been suggested by the purpose of ascertaining whether a 
link might exist between Eq. (5.1) and certain problems inherent in Yang-Mills 
theory and Differential Geometry (Ts), |l|l, 0. 

Precisely, we show that applying the reduction technique to the spherically 
symmetric version of Eq. (5.1), a connection can be established between so- 
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lutions of Eq. (5.1) and solutions of a reduced nonlinear ordinary differential 
equation, i.e. Eq. (5.15), which is invariant under conformal transformations. 
Equation ( |S.15[ ) comprises interesting cases, which come from n = 1, 3, 4, 6 
and 10 respectively. For instance, for n = 10, Eq. (5.15) can be interpreted as 
an extended (elliptic) form of the e quatio n governing the Thomas-Fermi model 
of an atom, while for n = 4, Eq. ( 5.15| ) turns out to be related to Eucl idean 
Yang-Mills equations via 't Hooft's ansatz [0. The properties of Eq. ( 5.15 ) 
can be reflected on those of Eq. (5T). This is illustrated in Sec. 6. 

Finally, Sec. 7 contains some concluding remarks, while in Appendixes A, 
B, C and D details of the calculation are reported. 



2 Group analysis 



To the aim of looking for the symmetry algebra of Eq. (1.2), we shall apply the 
standard procedure described in ||l^. In doing so, let us introduce the vector 
field 

F = ea, + r,a, + c9, + (2.1) 

where ^, rj, ^ and </> are functions of x, y, z, u, and = 4-, and so on. 



A local group of transformations is a symmetry group for Eq. (1.2), if and 
only if 

pr(2V[A]=0, (2.2) 

whenever A = for every generator V of G, where pr^'^^V is the second prolon- 
gation of V p| . 

Equation (2.2) give rise to a set of constraints in the form of partial dif- 
ferential equations which enable us to determine the coefficients ^, 7], C, and (j). 
This has been carried out in part using a computer, by means of the symbolic 
language of [0. However, to facilitate the understanding of the method to non 
specialist readers, in Appendix A we write pA^^V explicily. 

A general element of the symmetry algebra of Eq. (1.2) is 



V = {cix + C2y + C3)dx + {~C2X + ciy + C4)dy + {c5 + cez)dz + 2 r-(c6-Ci)9„, 

7-2 

(2.3) 

where ci, C2, cg are arbitrary constants. The expression ( |2.3| ) holds for any 
allowed value of 7 (7 7^ 0, 1), except for 7 = 
cases later. 



2. We shall come back to these 



From (2.3) we obtain the following independent generators of the symmetries 
ofEq. dTfr 



Vi = xdx + ydy 

V2 = ydx - xdy, 
^3 = 9,, 



(2.4) 
(2.5) 
(2.6) 
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Vi = dy, (2.7) 

V^^d,, (2.8) 

VQ^zd,^2^du. (2.9) 

These operators satisfy the commutation relations 

[■^^1, V2] = [Vi, V^] = [Vi, V^] = Vz] = [V2, V^] = 

= [V^, Va] = [1^3, ■^^S] = [V^3, "^^6] - ["t^4, V5] = [^4, 14] = 0, (2.10) 

[VuV^]^~V^, (2.11) 

[■t^i, ^4] --^4, (2.12) 

[V2, Vz] = 14, (2.13) 

[Vz, 14>] = "^Is. (2.14) 



We deduce that the operators (2^)-(2^) form the basis of a six-dimensional 



solvable Lie algebra C containing a three-dimensional ideal {V2, V3, V4} which 
is isomorphic to the algebra of _B2, the Euclidean group in the plane. The 
center of C is zero, while its derived algebra is {V2, V3, V^}. Furthermore, C 
admits the Levi decomposition L = Li \> £2, where Ci = {V^, V4, V5} and 
£2 = {Vi, V2, Vq}- This decomposition is not unique. We can easily write 
other decompositions of this kind) . The symbol l> stands for the operation of 
semidirect sum. (To keep the length of the paper reasonable, here we shall not 
explain the algebraic terminology used above. Anyway, the reader unfamiliar 
with the previous mathematical concepts, could address, for instance, to [|T8|). 

From the elements Vi and Vg we give the generator of the coordinate scale 
transformation, that is Vq = Vi + Ve xdx + ydy + zdz', V3, V4, V5 generate x, 
y and z-translations, V2 is a rotation symmetry operator, and Vi and Vg have 
the meaning of generators of two dilations together with a translation of 2:^^ 
along the =pM directions, respectively. 



As we have already mentioned, the symmetry generator (2.3) does not in- 
clude the value 7 = §• This is due to the fact that a special case emerges, jus t 



for 7 = |, in solving the equatio ns d et erm ining the coefficients involve d in (2.1). 



In this case, to the generators (2.4)-( p.9[ ) of the symmetries of Eq. (|l.2|) , one 



needs to add another operator, namely 

Vt = z'^dz- z{2 + ?,u)du, (2.15) 
which commutes with Vi, V2, V3, V4, while 

[V^.Vj]=Vt, [V^,Vj]^2V^. (2.16) 
A second special case to be dealt with separately is 7 = 2, in correspondence 



of which Eq. ( L2 ) becomes the linear wave equation 

k 
2 



yy-7^U,z=0. (2.17) 
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This equation is the continuous approximation of the lattice system arising from 



the harmonic potential 4>{rn) — jan.bnf''^ (see (1.3)). The symmetry algebra 



related to Eq. (2.17) is represented by the vector fields 



Gi = x{zd,^ + ydy) + \{x^ ~ y"^ + 2kz^)d^ - ^xud^, (2.18) 

G2^xd, + ydy + zd,, (2.19) 

G3 = xd, + 2kzd^, (2.20) 

G4 = ~xdy + yd^, (2.21) 

G5 = dy, (2.22) 

Ge = d,, (2.23) 

Gr = 9., (2.24) 

Gg = (2.25) 

Gg - fix, y, z)du, (2.26) 

where / is an arbitrary solution of the wave equation f^x + fyy — 2k f^^ — 0. 
The commutation rules fulfilled by Gi,...,G9 can be easily found. 

For brevity, the case 7 = f and 7 = 2 will not be further discussed. We 

remind only the reader interested in going deep into the case 7 = 2, that an 



exhaustive analysis of an equation of the type (2.17) is performed in ll 



3 Exact solutions from reduced equations 

The technique of symmetry reduction of a field equation amounts essentially to 
finding the invariants (symmetry variables) of a given subgroup of the symmetry 
group allowed by the equation under consideration jl^ . This method is usually 
preceded by the classification of the subalgebras of the symmetry algebra. The 
classification scheme is based on the adjoint action of the symmetry groupjjlS . 



Howe ver, because of the finite-dimensionality of the symmetry algebra ( 2.10 )- 
( ^.14| ), we shall adopt a more heuristic procedure. In other words, since for 
practical purposes generally one confines oneself to handle only low-dimensional 
symmetry subalgebras, which can be recognized in our case directly by inspec- 
tion, we do not need to employ the adjoint subgroup classification method. 

For each symmetry group generator V, we can obtain a basis set of invariants 
/(x, y, z, u) by solving the first order partial differential equation 

Vlix, y, z, u) = 0. (3.1) 



In the following, we shall build up the reduction of Eq. (1.2) by the one- 
dimensional subalgebras {Vi}, {V2}, {Vg}, and by the two-dimensional subal- 
gebras {Vi, V2}, {Fi, ^5}, {^2, Ve}. 
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3.1 Reductions by two— dimensional subalgebras 



The reduced equations associated with the two-dimensional subalgebras i) {Vi, V2}, 
ii) {Vi, V5} and Hi) {V2, Vq}, turn out to be nonlinear or dina ry differential 
equations which can be solved to give exact solutions to Eq. ( |l.2| ). 
Case i) 

A basis of invariants of the subgroup of the subalgebra {Vi, V2} is furnished 
by the system of equations Vil — 0, V2I = 0, which have to be satisfied simul- 
taneously. We obtain 



from which the reduced equation 



2fc7 



27 (7-l)(7-2) 



= 0, 



with A' = is derived. 



Putting 



into Eq. (3.3), we get 



A = 9- 



(3.2) 



(3.3) 



(3.4) 



7 



yielding 



7-1 V2 - 7 

(7-l)(7-2)2 

where c is a constant of integration. 
Equation (O) leads to the relation 



+c, 



(3.5) 



(3.6) 



1 + 66*'") ^ d0' = ±y/i{z - Zo) 



(3.7) 



zo being an arbitrary constant, and 
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c(7-l)(7-2)2 



(3.8) 



The integral at the left-hand side of (3.7) can be calculated resorting to the 
formula (see 111, p. 284) 



{l + bX'Y 



■2Fi{v, fj.; + -bX), 



(3.9) 
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where arg(l + bX) < tt, Re^i > and 2F1 denotes the Gauss hypergeometric 

I — A I — k 
function. Indeed, setting in ( p.Tl) X = 9'' and using ( |3.9D with IJ' — j and 

2 ' 



= i we find 



62 Fi 



I 2. 2 

2' 4' 4 



1; -6^?- =±^^(z-zo), 



(3.10) 



where 9 — 9(±y^{z — zq)) is exphcitly known whenever Eq. (3.10) is invertible. 
When this occurs, Eq. (|3.2|) provides an exact solution to Eq. (1.2). For 



example, for 7 = 4, fc = — 1 and zq — — |\/c, we have 



2F1 



-, 1; 2; -b9 



2(-i + vr+^ 



(3.11) 



with b = 

3c 



Then, combining together Eqs. (B.IO) and (3.2), we obtain the solution 



u = 4 



3c 1 
16 x'^ + 



to the equation 



^xx ~t~ ^yy ^" 



. 64 , 

1 z' 

9c 



1+4 



-4 



(3.12) 



0. 



(3.13) 



Case a) 

The subgroup of the subalgebra {Vi, V5} admits the basis of invariants 



y z 



(3.14) 



Then, the reduced equation of Eq. (1.2) coming from the symmetry variables rj 
and B{r]) takes the form 



B" 



B' 



7 



2-77 B 2-7 



with B' 

By introducing the change of variable 

B^ _ 2-j Q' 
B ^ 2 Q ' 



Eq. (3.15) can be cast into the linear differential equation 



rf{r]' + l)g" + 2ry rf + 



2-7 



Q' 



27 

(2-7)^ 



= 0, 
(3.15) 

(3.16) 
(3.17) 
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Resorting to the MATHEMATICA symbolic language, we find the general so- 
lution 



Q{v) = 



l+T]' 



,2 \ — 



ci sm 



2-7 



arctangry + C2 cos 



2-7 



arctangry 



(3.18) 

of Eq. (3.17), where ci and C2 are arbitrary constants. 

Thus, by virtue of ( 3.1^ ) and ( 3.16 ), an exact solution of Eq. (L2) is deter- 
mined, namely 

(3.19) 



u ~ J 



Case in) 

A set of basis invariants for the subgroup of the subalgebra {V2, Vg} is 



/i =T= (a;2+y2)^^ /2 = G(t)-z 1 + 



These give rise to the reduced equation 



2-7 



2-7 



of Eq. (|TJ), with G' = ^. 

Equation ( [3.21 ) can be written as 



via the transformation 



9" + -0' - 2fc 7 \ fl^-i 

T (2-7)^ 



2-"/ 



Thus, a particular solution of Eq. ( 3.22| ) is 

2k 

7 — 1 T'^ 



which yields (see (|3.20D and (|3.24|) ) 



M = 7 



2k t\ 2 

7 — 1 2; 



(3.20) 



(3.21) 



(3.22) 
(3.23) 

(3.24) 
(3.25) 



3.2 Reductions by one— dimensional subalgebras 

In opposition to what happens in the case of two-dimensional (symmetry) sub- 
algebras, the one-dimensional subalgebras provide reduced equations of (1.2) 
given by partial differential equations in two independent variables. Hence, in 
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order to arrive at reduced ordinary differential equations, we have to apply again 
the reduction procedure. This emerges in dealing with the one-dimensional sub- 
algebras i) {Vi}, a) {V2} and in) {Vq}. 
Case i) 

Associated with the symmetry vector field Vi , the basis of invariants 



X 1 ^ 

h ^ X, /2 = 77 = -, 13= A{ri, z) = - ( 1 + - 



(3.26) 



can be constructed. Then, substitution from the variables ( |3.26 ) into Eq. (1.2) 
yields the reduced equation 



7^2(772 + i)AA" + 277(772 + a)AA' + {a- 1)77^(772 + l)A'^ + {a - l)A^ 



k a-2 
2 a-1 



{a-3){^)+^ 



(3.27) 



where A' 



OA 



) Ax — and 



2 

,2-7, 



The symmetry algebra of Eq. ( 3.27 ) is expressed by the generators 

Ti = zdz - udu, 
T2 - (7?2 + 1)9, - 

T3 = a,. 



(3.28) 
(3.29) 
(3.30) 



The reduction by the vector field Ti leads to the invariants 

/l=7/, /2 = C/(77) = ZA(77, z), (3.31) 

which once replaced into Eq. ( ^.27 ) provide the ordinary differential equation 

772(772 + 1)UU" + 277(772 + a)UU' + {a- l)rf{rf + l)U'^ + {a- l)U'^ = " 2), 

(3.32) 

fied by putting 
U^e^. (3.33) 



with U' = This equation can be simplified by putting 



Indeed, we have 



772(772 + 1)9" + 277(772 + a)9' + a{a - 1)9 = -a(a - 2)9' 



(3.34) 



where the following 



2 a-1 



(3.35) 
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is a special solution. Now, taking account of ( ^.35| ), ( |3.33 ), ( |3.31 ) and (3.26), 
we find 

I- ] 

'7-1 2- 



u = 7 



7 



1 



(3.36) 



Furthermore, we observe that (3.32) allows the constant solution U = y f fiif : 
to which the solution (symmetric of (3.36) with respect to the exchange y — > x): 

'7- 1 x^^ ^ 



u — ^ 



7 z'' 



1 



(3.37) 



to Eq. (1.2) corresponds. 

Another interesting reduced equation can be derived from the symmetry 
generator T2. The related invariants are 



h =2, h = W{z) 



A{x, r])ri 



The introduction of (^^) into Eq. ( |3.27D gives 

W" + {a- 3)W' = 2fc^^^^W^3 



a- 2 



(3.38) 



(3.39) 



with W' 

For a = 3 (7 = I), Eq. ( 3.39 ) becomes 

W" = 12kW^, 

which affords the solution 



W 



1 1 



(3.40) 
(3.41) 



for k = —1, zq being an arbitrary constant. 

Another solution of Eq. ( 3.40D can be obtained carrying out a first integra- 



tion, which furnishes 



W' = 6kW^ + const. 



(3.42) 



A second integration provides a solution expressed in terms of elliptic functions 
(@, p.543). 

In general, by setting W — S,^ °o'(C); C = c^, Eq. (3.39) can be written as 



+ (4 - a)a, - 2fc^^-i^r^"+V3 = 0, 
that is an equation belonging to the class 



(3.43) 



(3.44) 
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(ai, 02 are constants) studied by Flower and by other authors, the references of 
them are quoted in p. 560). 

Case a) 

The symmetry subalgebra {V2} gives rise to the reduced equation 



T 



7-1 



where the basis of invariants 

Ii = T = \/ x'^ + y'^, I2 = u{t, z) 

has been used. 

A group analysis of Eq. ( ^3.45 ) provides the vector fields 

5*1 = zdz + Tdr 

S2 = ^Zd, + (J + U)du, 

S3 = d,. 



(3.45) 



(3.46) 



(3.47) 
(3.48) 
(3.49) 



We note that for 7 = |, in addition to Si, S2, S3 a further generator exists, i.e. 



'5'4= [ -+u] zdu~ -z^d, 



The reduced equation coming from (3.47) reads 



7 dp^ 



where 



u z 
V=l + -, p=-, u^u{p). 
7 r 



A particular solution of Eq. ( p.5l[ ) is 

2k \^ 



V 



which yields 



w = 7 



7-1. 
2k \ ^ /z\ 7^ 



7- 



(3.50) 

(3.51) 
(3.52) 

(3.53) 
(3.54) 



from ( 3.52 ), with r = \J + . 

On the other hand, the generator ( [3.46 ) leads to the reduced equation 



W" ^ -W' ^ XW^"^ , 

T 



(3.55) 
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where w{t) = r - v^^MV, W^' = ^ and A = ^j^- 



A solution of Eq. ( 3.55 ) is 

W = WoT-^, (3.56) 
with Wo = 7 ^ ^ . In correspondence of ( 3.56| ), the foUowing solution 



7 



(3.57) 



to Eq. (|L2D is found. 

Case in) 

By solving the equation Vg/ = 0, we get the set of basis invariants 



h 



X, h = y, h = W{x, y) ^ In 



1 / u 

— 1 + - 

z V 7 



7-2' 



Then, from (1.2) we have the reduced equation 



1 



yy 



2fc(7-l) w 



(3.58) 



(3.59) 



7-2^ ^ 7-2 
The group tecnique can be applied to obtain a further reduction of Eq. 



( p.59| ). In doing so, we find the symmetry generator 

iV = {-~c\x + C2y + cz)dx + {-C2X - ciy + Ci)dy + 2cidw- 

From ( [3.60 ) we deduce the independent vector fields 

Ni = -xd^ - y9j^ + 2(9^^, 
N2 = ydx - xdy, 
N3 = d,, 

N4 = dy, 

which fulfill the commutation relations 

[iVi,7V2]=0, [TVi, 7V3] =7V3, [TVi, iV4] = 



(3.60) 



(3.61) 
(3.62) 
(3.63) 
(3.64) 



(3.65) 
(3.66) 



It is noteworthy that Eqs. ( p.65| )-( |3.66| ) define a Lie algebra isomorphic to 
the algebra sm(2) of the similitude group in TZ^, governed by the rules pl| : 
[Xi, X2] = —X3, [Xi, X3] — X2, [Xi, X3] = 0, [X2, X4] = X2, [X3, X4] — X3, 
[Xi, X4] = 0. This can be checked by identifying A^i, N2, N3, N4 with —X4, 
—Xi, X2, X3, respectively. 
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Now let us examine the generator iVi , whose a set of symmetry variables is 



Ii=t^-, l2 = (T = 2lTix + W{x,y). 

X 



The reduced equation of (3.59) emerging from (3.67) reads 



1 



rtcFt 



27 2fc(7-l) 



which is solved by 



= In 



7-2 7-2 7-2 
fc7(a2 + 62) 1 



(3.67) 



(3.68) 



(3.69) 



7-1 {a + bty 

where a and b are consta nts o f i ntegr ation. 

Combining together ( 3.69 ), ( 3.67 ) and ( 3.58 ), we obtain the class of 00^- 
solutions 



w = 7 ■ 



fc7(a2 + 9) z2 



7 — 1 [ax + byY 



- 1 



(3.70) 



At this stage it is instructive to look for other nontrivial solutions of ( |1.2D . 
To this aim, let us set 

H'=(7-2)lnV' (3.71) 



into Eq. ( 3.59 ). Then, this equation takes the form 

'4^xx ~t~ '4^yy ^ 2/c 



'(7-2)2 



(3.72) 



By way of example, we seek a solution of the type • = '0 (0, where ^ ~ x — vy 
and w is a (real) constant. In such a manner Eq. ( 3.72| ) transforms into the 
(nonlinear) ordinary differential equation 



l + i;2(7-2) 



from which 



-0? = aip^ + c. 



(3.73) 
(3.74) 



where c is a constant of integration and a = 7(7-2)^ ' c — Q Eq. ( 3.74 ) 
provides 



V' = ± 



a2 



2-7 



(3.75) 



^0 being a constant of integration. 

On the other hand, if c 7^ we obtain 



^1 + b^'^ 



- &) 



(3.76) 
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with 6= 

The integral at the left-hand side of (3.76) can be evaluated by Eq. (3.9). 
Indeed, putting in (3.76) X = tjj^ , from Eq. ([3.9[) we find 



7^2Fi ( -, -; 



(3.77) 



with LL = - and ly = h 

^7 2 



With the help of ( 3.58| ) and ( 3.71 ) , we get the exact solution 

2 

u — "/[zi-^ip — 1) 



(3.78) 



t o Eq . (1.2), where tp 



ip{x — vy) is explicitly known if one can invert Eq. 



Equation ( 3.77 ) covers some cases related to potentials of physical in tere st, 
arising for a) 7 = 3, b) 7 = —2, c) 7 = |, d) 7 = |, and e) 7 = — 1 (see (1.3)). 

The first two choices correspond, respectively, to the Fermi-Pasta-Ulam 
potential |^ and to a potential whose nonhnear part, 1/ (l — f ) , finds applica- 
tions in the treatment of the scattering states in conformally invariant Quantum 
Mechanics In cylindric coordinates, case c) leads to an equation involved 
i n the Thomas-Fermi model of an atom (||l^ p. 116), while in the case d) Eq. 
( ^.72| ) appears in the theory of the white dwarf stars ((l^, p. 113). Finally, for 
-1 (case e)) the potential (L3) takes the form (/> ^ (\ — u)^^ — {l + u\ which 
mimics a special case of the Killingbeck potential V — + Br + Cr^ |Q. We 
recall that a Coulomb potential perturbed by a second degree polynomial can 
be used to study the ground state properties of a hydrogen atom. 
For brevity, below we handle in detail only the cases a) and b). 

Case a) 

For 7 = 3, Eq. ( ^.74 ) can be integrated straightforwardly. Hence, first we 
shall evaluate tp explicitly from ( 3.74 ). 

Subsequently, we shall compare the expression of ?/; determined in such a 



7 



way with that coming from (3.77). In doing so, let us reshape Eq. (3.74) into 



93, 



(3.79) 



after the rescaling i] 



^ ^, with c : 



3(l + D^) 



-3(1+^93 and k 
Equation ( [3.79| ) is a special version of the equation 



(3.80) 



which is satisfied by the Weierstrass elliptic function p(ry; 52, 53), where (72 and 
(73 are the invariants of p (|2j] p. 629). 
For 93 — we have 

^(77) = Pirn 0, 0) = 4' (3-81) 
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that is the first term of the series representation of p{ri; g2, gz) ([Q p. 635). 
Then, (|3J|) gives 



u = 3 ^ - 1 



77- 



For 7^ 0, Eq. (3.8C) yields 



namely 

Thus, we infer that 



■0 = p{m 0, 53)- 

U = 3 [z2p(,7; 0, 53) - 1] 



(3.82) 

(3.83) 

(3.84) 
(3.85) 



from (|3.78| ). 

Combining ( 3.84 ) together with ( [3.77] ), it can be shown (see Appendix B) 
that the remarkable formula 



1 1 4 



p-' = ^2Fi -, -; -; -4^3 - —B -, - 



2' 3 3 



41/3 



1 1 



3' 6 



(3.86) 



holds, where p ^ stands for the inverse Weierstrass function (92 = 0), = — ^ 

and B(p, q) is the beta function (Euler's integral of the first kind) ([191 p. 
948). The relation (3.86) estabilishes a link between the Gauss hypergeometric 
function and the inverse Weierstrass function (92 = 0). 
Case h) 

If 7 = -2, Eq. ( ^ ) furnishes 



where b — 



8k 



3c(l+i,^)' 



1 1 1 



-btlj- 



(see [!§, p. 1042). 
Consequently, 



u = -2 



(3.87) 



(3.88) 



(3.89) 



from (|3J8|) . 
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4 Instanton— like approximate solution 



Equation (3.59) can be exploited as well to give approximate solutions to Eq. 
( L^) of the instanton-like type. To this purpose, it is convenient to write Eq. 
( 3.59| ) in the form 



dr,djfW - 



27-2 ' 



(4.1) 



where rj — x + iy, 1] = x ~ iy. 

For 7 ^ CX3, Eq. (iA) reproduces the Liouville equation 



dr,drfW = -e 



w 



(4.2) 



which is conformally invariant, in the sense that if W{ri,ri) is a solution of Eq. 
(^), the function W = W{'ij, 77) given by 

W^W-lnf'j' (4.3) 

is also a solution, where r/ = /{rj), 77 = /(jy), and /' = ffjif]), / = f-ff{fi). 

We notice that the Liouvil le e quation emerges as a symmetry reduction of 
the continuous Toda equation (LI) |lj] corresponding t o th e generator V6(oo) = 
zdz + 2du, coming from ( ^.9| ) for j ^ 00. Equation (L2) affords the general 
solution _ 

[1 + fiv)fivW 

for k = —1, where /(r;) denotes an arbitrary holomorphic function. This ex- 
pression can be exploited to get (exact) instanton solutions to Eq. (1.1) ]l^ , 
[14|. Conversely, we are not able to obtain exact instanton solutions to Eq. 
( [4.1| ). However, we can find (heuristically) an approximate instanton-like solu- 
tion, which arises from the assumption 







(4.5) 



where Wq satisfies (L4), e = is such that |e| ^ 1, and </) is a function of 77, 

77 to be determined. 

Substitution from (L5) into (|4.lD gives 



W, 



Orjrj 



(4.6) 
(4.7) 
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for large values of 7. At this point we choose k = —1, / = 77, / = r; in (4.4). 
Then, Wo becomes the instanton configuration 



(1 + 7777)^ ' 



With the help of (4^), Eq. (4/7) can be written as 



g„ + -g^ + — c 
r 



8.9 



(l + r2)2 

in polar coordinates (a; = rcosO, y — rsinO), where 

g{r,e) ^ (t){r,e) + 2r^ + 1. 



(4.8) 



(4.9) 



(4.10) 



We are interested in solutions of (|4.9| ) of the type g = g{r, 0). Therefore, by 

/ n1/2 

introducing the change of variable r = ( j , Eq. (|4.9| ) takes the form 



(1 - t^)9tt - 2tgt + 2g 



{1-ty 



(4.11) 



It is noteworthy that, in general, this equation can be solved exactly. In fact, its 
homogeneous part coincides with a special case of the hypergeometric equation 
defining the Legendre functions (Q, p. 331). Precisely, the homogeneous part 
of (4.11) allows the independent integrals 



gi{t) ^ Pi{t) = t, 
92{t) = Qi{t) 



1 , t + l 
-t\n 1, 

2 1 



(4.12) 
(4.13) 



where Pi and Qi are the Legendre polynomial and the Legendre function of the 
second kind, respectively, corresponding to tt, = 1. 
In terms of , we have 

31 



1 



92 = 7; 



+ 1 

„2 



2 r2 + 1 



Inr — 1 H gi. 



Furthermore, a particular solution to Eq. (4.11) is 

o 7-2-1 



5o = 2<i41n(l + r- 
where 



- [2dilog(l 



Inr 



dilog(x) 



Int 



-dt 



(4.14) 
(4.15) 

(4.16) 
(4.17) 
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is the dilogaritm function ([Q p. 1004). 

Now coming back to the original function ip (see (|4.1C )), we are led to the 
following instanton-like approximate solution to Eq. (4.1) (see (4.5), (4.8)): 



W = In 



(l + r2)2 



- 1 



e < ci- 



C2 



1 - 1 



2r2 + 1 



In r' 



-1 



+ 81n(l + r2) - 10 + 2— f2dilog(l + r^) + Inr^l - iL (4.18) 

+ 1 '■ ^ } 

ci, C2 being arbitrary constants. This formula is meaningful for small e (large 
7). 

In order to confer an instanton character to the solution (4.1J), we need to 
choose C2 = —4. Thus, the function W becomes singularity free, precisely: 



W = In 



(l + r2)2 



[ci + 4dilog(l + r2)] +81n(l + r2)-7|. (4.19) 



For small values of e, the function ( 4.19|) can be used as a good approximation 
of the regular solution Wi (r) of Eq. ( 3.59| ) expressed in terms of r. In Appendix 
C we report two Tables, la and lb, containing values of W{r) and Wi{r) for 
some r with 7 fixed. The initial conditions have been chosen in such a way that 
ci = — 7 and C2 = —4. From the Tables one argues that W{r) and Wi{r) are 
quite closed already for 7 = 6 (i.e. e — 0.25). The quantity |VK(r) — VF/(r)| 
tends to zero as 7 increases. 



5 Extension of Eq. (ri.2|) to (n + l)-dimensions 

Let us consider the generalized version 



(92^+...+a2J u = k 



1+^ 

7 



7-1' 



(5.1) 



of Eq. (1.2), where M = M(a;i, a;„, z), with (xi, a;„) S 7?.". 

Assuming that u — u{r, z), where r = •^a;2 + • • • + a:2 , Eq. (5J) can be 
written as 



n- 1 



1 + ^ 

7 



7-1 



(5.2) 



The application of the symmetry approach to Eq. (5.2) provides the symmetry 
generator 



G = kirdr 



7-2 , 
fci H 7: — k2 + fcs 



dz + fc2(u + 7)(9„, 



(5.3) 



with ki, k2, fcs arbitrary constants. 
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From ( p. 3D we get the operators 

Gi = rdr + zdz, 
7-2 

G2 = —^zdz + (7 + u)du 
G3 = d,, 
which satisfy the commutation relations 

[Gi,G2]=0, [Gi,G3]=— G3, [G2,G3] 



7-2 



G, 



(5.4) 
(5.5) 
(5.6) 

(5.7) 



The symmetry operator G2 is of particular interest. In fact, it gives rise to the 
invariants r' — r, W'{r') = W{r), with 



W{r) = z T-a (-y + u), 
which leads to the ordinary differential equation 



Wrr 

Equation (^.9[) becomes 



r " 77-2(^-2)2'' 



Iprr H tpr = Qlp'^^^ , 

r 



(5.8) 



(5.9) 



(5.10) 



with g — '^^^_2y , by setting W — "/tp. The link beetween u and ip is given by 
(see (0)) 



u = "f (^z-1-^4' ~ 1 
Now we look for a transformation of the type 

r' = -, = /(r')^(r'), 



which leaves Equation (5.10) invariant, in the sense that 



n-1 



r 

A straightforward calculation shows that this can be accomplished if 



r/ f\ /n~2 



2n 
n- 2' 



In this case Eq. ( 5.10 ) takes the form 



n-1 r, + 2 

Wrr ^ Wr = SoV'""^ : 

r 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 
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with 



k 



go-|(n2-4) (n^2). 



(5.16) 

Furthermore, it is easy to see that if ^p{f) is a sohition of Eq. ( 5.15| ), i.e. 

n — 1 



Iprr 



(5.17) 



then ^{r) = \^~^p{r), where r = \r and A is a real parameter, is also a solution. 

To summarize, we have 
Proposition. Equation (5.15) is invariant under the conformal transformation 



r' = — , ^jJ(r) — r'" 2 ^p(r'). 
Ar 



(5.18) 



The positive (spherically symmetric) regular solution to Eq. (5.K) is given 



by 



-0(r) 



(l + r2)^ 



(5.19) 



with a = (^—-^^^ , where k and n are to be chosen in such a way that a > 0. 

We notice that Eq. ( 5.15 ) comprises remarkable special cases, which are 
listed in the following Table: 



Table II 



Case 


n 


7 


Equation (5.15|) 


Regular solution il>{r) {k = —1) 


a 


1 


-2 


Iprr = -^klp-^ 


(_^)-^/* (1+^2)1/2 


b 


3 


6 




(_^fc)'/"(l + r2)-i/2 


c 


4 


4 


Ar + f V'r = ^kA 


(-ffc)^/^(l + r2)-i 


d 


6 


3 


Ar + ^A= 4^1/-^ 


-6fc(l + r2)-2 


e 


10 


5/2 




^(l + r2)-4 



The solutions of the original Equation (5.9), related to cases a), 
from 

2n 



e), arise 
(5.20) 



(see (pTJ)). 

It is noteworthy that Eq. (5.15) appears in the context of Differential Geom- 
etry ||l^ , jl^ , in particular it is related to the socalled Yamabe problem, which 
consists essentially in establishing when a Riemannian metric can be changed 
by a conformal (lenght) factor to have constant scalar curvature [E5|. 
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Moreover, some special cases of Eq. ( 5.15| ) play a basic role in certain 



branches of physics. Precisely, the equation of case b) of the Table II con- 
stitutes a static and spherically symmetric version of a nonlinear wave equation 
discussed by Rosen some years ago Equation b), which is of the Emden 

type, finds application in Astrophysics. 

Equation e) can be interpreted as an extended (elliptic) conformal invariant 
version of the equation 

^rr + -*r = C^^^^ (5.21) 

r 

governing the Thomas-Fermi model of an atom, where C > is a certain 
constant and i'(r) is the potential field originated by Z — 1 electrons, acting on 
the Z^"^ one ([...], p. 125). 

Equation a) corresponds to a potential whose nonlinear part, 1/(1 — ^)^, 
mimics the inverse square potential appearing in the treatment of the scattering 
states in conformally invariant Quantum Mechanics Finally, Eq. d) is 

associated with the Fermi-Pasta-Ulam potential. 

To conclude the group analysis of Eq. (|]^), wc mention another nontrivial 
reduction, which can be written down starting from the generator Gq = Gi + 
G2 + G3. It reads 

X^Wxx + {n- 1){W - XWx) = ^^x^W^-^WxX + 

Ih^my^wz_^^2XWx + W), (5.22) 



where X = undW ^W{X) ^ ^ are two invariants. Equation (|^ 

will not be discussed here. 



6 Relation to the Yang— Mills theory 



Equation c) of Table II can be considered as a reduction form of the Yang-Mills 
equations (see later). This is a well-known result ||l^, [Q. What is new, is the 
link between special solutions of the Yang-Mills equations and special solutions 
of the field equation (5.1) for n = 4. To this aim, let us start from Eq. (5.10), 
which can be written as 



n-2/n-2 \^ n - 2 

Clrr H a fir - 2- 



a-1 \n-2 
via the transformations 



(a-iy 



a 



71-2 



1 ^ 



n + Q°'=0, (6.1) 



(6.2) 



with 11 = 2a ' ^^^"^ Q = 7 — 1. We notice that Eq. (6.1) is the same 

as Eq. (1.6) contained in |l^]. Therefore, all the results achieved about this 
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equation can be transferred to the field equation (5.1) through the formula 
( ^.11 ), namely 

' a — 1 



M = (a + 1) 



2a r 



(6.3) 



where and r are expressed by (^]^). 

In general, Eq. (6.1) may be investigated using phase-space variables. 

Here we restrict ourselves to the case a = in correspondence of which 
Eq. (6.1) reads 



n- 2 



n + n^ ^ 0. 



Equation ( |6.4[ ) can be solved exactly. Indeed, we easily find 

2 



n - 2 



H r2"-2 



where c is a constant of integration. 

By choosing for example c = 0, a simple calculation gives 



n{r) 



\J n{n — 2) • r 



A2 



or, in the variable ip{r) (see (6.2)): 



tp{r) = 



4A 



1 



2(n + 2) r2 + A2_ 
This function represents a regular solution to the equation 



71 — 1 
1p„ H T/^r 



71^ — 4 , >.+2 



8 



=0, 



(6.4) 



(6.5) 



(6.6) 



(6.7) 



(6.8) 



which comes from Eq. ( |5.10| ) for a = 7 — 1 = On the other hand, replacing 
(3.6) into (|6.3|) provides the exact regular solution 



271 

71-2 



4A 



2(71+2) z2 + A2 



(6.9) 



to the field equation ( pTll ) (^7 

Equation (3.4) admits also singular solutions The link between this 

kind of solutions and the corresponding ones of Eq. (5.1) can be obtained via 
(pT^). At this stage, we point out that for rt = 4, a connection can be established 
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between solutions of Eq. (^J) and solutions of the Euclidean Yang-Mills (YM) 
equations. To pursue this goal, let us introduce some preliminaries. 
The pure YM equations (i.e., in absence of matter fields) are 



(6.10) 



(fi v = 0, 1, 2, 3), where the summation convention is understood. 

The fields (x) take values on the Lie algebra G oi a compact semisimple 
Lie group G (the gauge group). By choosing SU(2) as the gauge group, we can 
write 

(6.11) 



(a = 1, 2, 3), where are the Pauli matrices and A'^ = iTr{Ap(a;)(T°}. 
On the othe hand, the tensor i^^,^ associated with A^i(a;), defined by 



F^Ax) 



d^Afj^ + [A^, A^] 



is expressed by 



F^.,{x) 



1 



(6.12) 
(6.13) 



^it\iF-Ax)^iTY{F^,{x)a-}. 

Now, inserting the 't Hooft ansatz |17 



-11atLydy\n^j{x) (6.14) 

into Eq. ( |6.f 0| ), where the tensor rya/^iywill be specified later (see Appendix D), 
we obtain the equation 

VV + -ft'V'^ = 0, (6.f5) 

K bei ng a constant. By writing in spherical coordinates and putt ing K — | , 
Eq. (|]f|) becomes just Eq. (|!|) for n = 4. Then, by virtue of dJ), 
(x£) and (O), from ( 6.14 ) we find the solution to the Yang-Mills equations 

= -Vat^pdy In 



i(l 



where u is the regular solution 



4A 



A2 



to the nonlinear field equation (see (5.2)) 

3 



(6.16) 



(6.17) 



(6.f8) 



with r = y/x-f + X2+ xl+ x\. 

For completeness, in Appendix D a detailed proof of the correspondence 
between Eq. (6.f 5) and (6.12) is presented. 
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7 Conclusion 



Using a group-theoretical approach, wc have investigated the nonHnear field 



equation (1.2), arising from a la ttice of the binomial type as a continuous approx- 
imation, and its extension to (n + l)-dimensions. For 7 ~> 00, both equa- 
tions ta ke the f orm o f Toda field equati ons . We have determined the symmetry 
algebra ( 2.10| )~( ^.14 ) admitted by Eq. (2^). This algebra is finite-dimensional, 
while that allowed by the Toda field equation (LI) is infinite-dimensional. This 
result constitutes a first discrimination between the two equations. Another 
important difference is that Eq. (LI) enjoys the confor mal invariance property. 
On the contrary, this property is not valid for Eq. (1.2), but it is regained 
through Eq. (5.15) (n ^ 2). In other words, the conformal invariance property 
holds again within the generalized equation (5.1). This feature enables us to 
build up conformal versions of physically interesting models, listed in Table II. 
These conformal models afford both regular and singular solutions reflecting on 
the solutions of Eq. (5.1) via the transformation ( 5.20[ ), which comes from the 
reduction procedure applied to Eq. (5.2) (the spherically symmetric form of Eq. 
(^n|)). The introduction of the generalized equation ( |5.l| ) leads to a scenario 
where several nonlinear field equations appearing in different physical and dif- 
ferential geometrical contexts can be dealt with in a unifying manner. This is 
the spirit of the present work, which has essentially a speculative character. 



Appendix A: explicit form of pA'^^V 

The second prolongation at the left-hand side of Eq. (|2.2|) is given by 



where 





prWy ^ 












dU-cy 








syy^S- 

^ OUy 


y UUyz 






(A.l) 




)v = v 


IT d 
mix 


dUy 


duz 


(A.2) 


r = 




- VUy - (Uz) 


+ ^Uxx 4 


TJUxy H~ C'^xz 


(A.3) 






- T]Uy ~ (Uz) 


+ i,Uxy 4 


TjUyy ~h C,Uy2, 


(A.4) 


r = DM 




- rjUy - (uz) 


+ ^Uxz 4 


rjUyz + (Uzz 


(A.5) 




- - 


rjUy - C,Uz) -1 




1]'^xxy ~l" Q'^xxz 


(A.6) 


= D,Dy{ 


(j) - S^Ux 


- 'quy - Cuz) 


^'^xxy 


-f- fJUxyy ~t~ C'^xyz 


(A.7) 


.-^ = D^D.i 


(j> - ^Ux 


- -qUy - (uz) 




-|- TjUxyz "4^ Q'^XZZ 


(A.8) 


^yy = Dli^ 


- S,Ux - 


i]Uy - Cuz) ^ 


^'^xyy ^" 


- TjUyyy + Q^yyz 


(A.9) 


^y'' = DyD,{ 


(j) - i^ux 


- 'qUy - Cuz) 


^'^xyz 


-|- TjUyy2, ~h C'^yzz 


(A.IO) 


0- = Dl{<j> 


- S.Ux - 


rjUy ~ (uz) ~\ 




TjUyzz C'^ZZZ 


(A.ll) 
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and Dx, Dy, Dz denote the total derivative operators 



d 



d 



d 



d 



ox ou oy ou oz au 

Appendix B: derivation of Eq. (^TS^) 



(A.12) 



Let us start from ( ^79| ), which can be written as 



Vl + 4cr3 

1 

via the substitution il) = —g^a. Then 

da' 1 r x-^/^x 1 r' x-^/Hx 



where X = and ry is a constant. 

Now, we recall that ([fol, p. 285, formula 3): 



(B.l) 



*.9|(^/-'7oo), (B.2) 



[i + bxy 



(B.3) 



with |arg6| < tt, i?ej^ > i?e/i > 0, where B{x, y) is the beta function (Euler's 
integral of the first kind) defined by (ll^, p. 948) 



B{x, y) = / e-\l-t)y-^dt. 
Jo 



(B.4) 



Taking account of pjj ) and ( |B.4| ) with 6 = 4, ^ = i and v = ^, Eq. ( |B.2D 
implies 



da' 



Vl + 4:a 



i3 



3 V3 6 



1 1^ 4^ 

2' 3' 3' 



-Aa' 



(B.5) 



Since a — p{ig§ {Vcc—il)'} 0, —1) (see ( 3.83 ) and ( 3.84 )), we obtain ig§ (?/oo— ??) 
p~^{a) and, therefore, formula (3.86|). 



With the help of the property B{x, y) 



be calculated explicitly. We have 
'1 1 



B 



3 6 



r(i)r(i] 



T{x+y) ' 



r(|! 



the quantity B \ 



^.413 



can 



(B.6) 



where the duplication formula for the gamma function, T{2z) = -^=2?-^ 2 
r(z)r (z+ i), has been applied to (z = i) ([||, p. 255). 
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Appendix C: some values of W{r) and Wi{r) for 
fixed 7 



Table la 



7 


1^/(0.1) 


T4^(0.1) 


1^/(0.3) 


W{0.3) 


1^/(0.5) 


W{0.5) 


4 


1.100 


0.763 


0.597 


0.406 


0.163 


0.069 


6 


1.148 


1.048 


0.729 


0.678 


0.369 


0.345 


8 


1.164 


1.117 


0.773 


0.750 


0.438 


0.427 


10 


1.172 


1.145 


0.795 


0.782 


0.572 


0.466 


12 


1.177 


1.159 


0.809 


0.800 


0.493 


0.489 


oo 


1.196 


1.196 


0.862 


0.862 


0.575 


0.575 



Table lb 
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Wi{l) 


W{1) 


Wi{2) 


W{2) 


Wi{3) 


1^(3) 


4 


-0.727 


-0.727 


-2.041 


-2.438 


-3.030 


-5.840 


6 


-0.364 


-0.364 


-1.426 


-1.487 


-2.208 


-2.418 


8 


-0.242 


-0.242 


-1.221 


-1.245 


-1.934 


-2.012 


10 


-0.182 


-0.182 


-1.118 


-1.131 


-1.797 


-1.837 


12 


-0.145 


-0.145 


-1.057 


-1.065 


-1.715 


-1.740 


oo 








-0.811 


-0.811 


-1.386 


-1.386 



Appendix D: solutions of Eq. (|t).lU|) via solutions 
ofEq. ( ^TTBI) 

The tensor ?7a^^ appearing at the right-hand side of ( |6.14 ) is defined as 



Vaij = ^ai], VavO = -VaOiy = Sau, VaOQ = 0, (D.l) 

where Caij denotes the Ricci tensor. Other properties of the symbol rj are re- 
ported in ijl^ . It is also convenient to use the matrices 

<7,j.u = ]^Vat,u(y'' (D.2) 

(here a, i, j, . . . — 1, 2, 3, and ii, v, \, . . . — 1, 2, 3). 

In order to write F^i, in ter ms o f cr^^,_we shall evaluate the commutator 



[cr^i/, Opx]. Taking into account ( D.l ) and ( |D.2| ), we find 



[Cpi/, CpA] = -r]at2iy'ribpX<^abcCr'^ , (D.3) 
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where the commutation rule [cr", a^] = 2i(T^ has been employed. On the other 
hand, we have 



Substituting (^) into (|DJ) yields 



(D.4) 



(D.5) 



with the help of (D.2) 



Now, the YM field ( |6.1lD can be written as (see ( |6.14| ) and (^J)) 

Ap = iapydy\iiil>. (D.6) 

Hence, the tensor Fp^, (see ( |6.12 )) takes the form 

Fp^ = icrj,p9^(9pln^) - iap^d^idplinp) ~ i{dplinp)idx\iii>)[app, a^x] = 

i{App - ApAp)a^p - i{A^p - A^Ap)(Tpp - i{Apfap^, (D.7) 

where the shorthand notation 

A, = 9,lnV^, A^p = d.dplniP, Al ^ dlln^j, (A.f = (9, In 7/>)2, ... (D.8) 

is used. 



K eepin g in mind ( |D.7| ) and (D^), the quantities dpFp^ and [Ap, Fp^] (see 
Eq. (6.12)) can be elaborated as follows: 



dpFp^ = i{Appp - AlAp - ApApp)(jyp 
-i{Apyp - ApyAp - AyApp)app - ^iApAppOpi. 



(D.9) 



{Ap, Fpy\ = iAp[crpp, Fpi,] 
-{ApApx - ApApAx)[app, ct^a] 

-(ApA^A - ApA^AA)[crpp, CTpA] 

+Ap(Aa)^[o-pp, ap^]. 



(D.IO) 



Then, putting (D.9) and ( D.10[ ) into the YM equations (3.10) and exploiting 
the commutation rule (D.5), after some manipulations we obtain 



a^p[Appp - A2 Ap + 2AaApa - ApA^ - 2Ap(Aa)2] 
+cjpp{-Ap^p + 2ApA^Ap) = 0. 



(D.ll) 



Since the tensor multiplyin g a^p is symmetric under the exchange ii p, while 
app is antisymmetric, Eq. ( |D.ll ) becomes 



Appp + 2ApA^p - 2Ap[A^ + (A^)^] = 0, 



(D.12) 
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which can be reduced to the form 



do hi 



InV/ 



(D.13) 



Equation ( D.13 ) tells us that the function in the brac ket is independent from 



Xp for any value of p = 0, 1, 2, 3. Consequently, Eq. ( p.lSD arises (9^ = V ), 
where X is a constant of integration. 
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